Loop quantum cosmology of the Bianchi I model: 
complete quantization 



M Martin-Benito^ L J Garay^ '^, G A Mena Marugan^, 
E. Wilson-Ewing^ 

^ MPI fur Gravitational Physics, Albert Einstein Institute, Am Miihlenberg 1, D-14476 
Potsdam, Germany. 

^ Instituto de Estructura de la Materia, lEM-CSIC, Serrano 121, 28006 Madrid, Spain. 
^Departamento de Ffsica Teorica II, Universidad Complutense de Madrid, 28040 Madrid, 
Spain. 

■^Centre de Physique Theorique, CNRS-Luminy Case 907, 13288 MarseiUe Cedex 09, France. 

E-mail: mercedes@aei.mpg.de, luis.garay@fis.ucm.es, menaSiem. cf mac . csic . es , 
wilson-ewingOcpt . univ-mrs . f r 

Abstract. We complete the canonical quantization of the vacuum Bianchi I model within the 
improved dynamics scheme of loop quantum cosmology, characterizing the Hilbert structure 
of the physical states and providing a complete set of observables acting on them. In order to 
achieve this task, it has been essential to determine the structure of the separable superselection 
sectors that arise owing to the polymeric quantization, and to prove that the initial value prob- 
lem obtained when regarding the Hamiltonian constraint as an evolution equation, interpreting 
the volume as the evolution parameter, is well-posed. 



1. Introduction 

Loop quantum cosmology (LQC) [H [21 [3] adapts the techniques of loop quantum gravity [3] 
in the quantization of models with high degree of symmetry, such us homogeneous models. 
Remarkably, the quantization of (homogeneous and isotropic) Friedman-Lemaitre-Robertson- 
Walker models within (the improved dynamics of) LQC succeeds in solving the singularity 
problem: the classical big bang turns out to be replaced by a quantum bounce happening at 
Planck scales, and no observable diverges in the quantum theory, as shown for the first time in 
the seminal work of reference [5]. 

In order to extent this quantization to (homogeneous but anisotropic) Bianchi models, we 
have focused our attention in the simplest case: the vacuum Bianchi I model. Our interest in 
this model also comes from the necessity of including inhomogeneities in LQC. Actually, with 
that aim we have analyzed from the LQC perspective the simplest inhomogeneous cosmologies: 
the linearly polarized Gowdy model with three-torus topology [6, 7J. This model can be regarded 
as a homogeneous Bianchi I background with three-torus topology filled with linearly polarized 
gravitational waves traveling in a single direction. The above quantization of the Gowdy model 
employs the polymeric quantization of the Bianchi I model when representing the homogeneous 
background and, therefore, first we need to have under control the Bianchi I model itself. 



The quantization of the Bianchi I model in LQC is subject to several ambiguities. One of 
them concerns the representation of the curvature tensor of the connection in the quantum 
theory. Different definitions of this object lead to different schemes of quantization, being 
the most satisfactory one the improved dynamics put forward in Ref. |8j, where although the 
Hamiltonian constraint was constructed, no further analysis of the physical solutions was carried 
out. 

In this note we further study this quantization of the Bianchi I model. It summarizes some 
results already presented in E] • We will start by adopting a different factor ordering of that 
of [5] when symmetrizing the Hamiltonian constraint operator. With our factor ordering the 
operator is well-defined in the octant of given by the positive eigenvalues of the operators 
representing the coefficients of the densitized triad. In this way our model displays sectors 
of superselection that are simpler than those of [S]. We will analyze the structure of the 
superselection sectors (aspect not studied in [8]) and see that they are separable. 

The Hamiltonian constraint provides a difference equation in the volume. Then, it seems 
natural to interpret it as an internal time and the constraint as an evolution equation with 
respect to it. We will see that this interpretation is valid inasmuch as the corresponding initial 
value problem is well-possed: physical solutions are completely determined by a countable set 
of data given in the initial section of the volume (which displays a non-vanishing minimum 
value). The previous result allows us to identify the physical Hilbert space with the Hilbert 
space of initial data, whose inner product is determined by imposing reality conditions in a 
(over) complete set of physical observables. 



2. Quantization of the model 

We consider the vacuum Bianchi I model with three-torus spatial topology. Then, we use global 
coordinates {t, 6, a, 6}, with 6,a,5 £ . In order to prepare the model for its loop quantization, 
we describe it in the Ashtekar-Barbero formalism Using a diagonal gauge, the nontrivial 
components of the densitized triad are pi/Air'^, with i = 6,a, 6, and Ci/ln are those of the su{2) 
connection. They satisfy {ci,pj} = SnG'jSij, where 7 is the Immirzi parameter and G is the 
Newton constant. The spacetime metric in these variables reads 

ds^=-dt^+\-^('4+'4+'4). (1) 

4vr2 \pj pI pj J 
The phase space is constrained by the Hamiltonian constraint 

2 c^pBC^Pa + c^pgc^ps + c^p^c^ps 

Cbi = -:^ V = 0- (2) 



In this expression, V = \J \pePaP5 \ is the volume of the compact spatial sections. 

In order to represent the phase space in the quantum theory we choose the kinematical 
Hilbert space of the Bianchi I model constructed in LQC (see e.g. [B]), that we call ^kin- We 
recall that, on ^kim the operators pi have a discrete spectrum equal to the real line. The 
corresponding eigenstates, \pe:Pa,Ps): form an orthonormal basis (in the discrete norm) of "Hkin- 
Owing to this discreteness, there is no well-defined operator representing the connection, but 
rather its holonomies. They are computed along straight edges in the fiducial directions. The 
so-called improved dynamics prescription states that, when writing the curvature tensor in terms 
of holonomies, we have to evaluate them along edges with a certain minimum dynamical (state 
dependent) length /ij. We use the specific improved dynamics prescription put forward in [8]: 
the elementary operators which represent the matrix elements of the holonomies, called A/^ij, 
produce all a constant shift in volume V. In order to simplify the analysis, it is convenient 
to relabel the basis states in the form \v,Xa,Xs), where \v\ is proportional to V such that the 



operators A/±^^ cause a shift on it equal to ±1, and the parameters Aj are all equally defined in 
terms of the corresponding parameters pi, and verify that v = 2XgXa-Xs- 

Out of the basic operators pi, Mfi^ we represent the Hamiltonian constraint as an operator 
C. We choose a very suitable symmetric factor ordering such that C decouples the states of 
T^kin with support in v = 0, namely, the states with vanishing volume [6l [7j. Moreover, our 
operator C does not relate states with different orientation of any of the eigenvalues of the 
operators pi. Owing to this property we can then restrict the domain of definition of C to e.g. 
the space spanned by the states \v > 0, Ao- > 0, Xs > 0). We call the resulting Hilbert space Ti^^^. 
Remarkably, for this restriction we do not need to impose any particular boundary condition or 
appeal to any parity symmetry. Note that both 'Hkin and V.^^^ are non-separable Hilbert spaces, 
feature not desirable for a physical theory. This problem is overcome by the own action of the 
Hamiltonian constraint operator. Indeed, C defined on ^^i^ turns out to leave invariant some 
separable Hilbert subspaces that provide superselection sectors. In concrete, they are spanned 
by the states \v = e + 4n, A^- = A^w^, A^ = A^a;^). Here e £ (0,4] and n G N, and therefore v 
takes support in semilattices of constant step equal to 4 starting in a minimum non-vanishing 
value £. In addition. A* (a = a, 6) is some positive real number and cug runs over the subset of 
given by 

where mk,nk,Pk G N, and z G Z when e > 2, while z = otherwise. One can check that in 
deed this set is dense in M"*" and countable [6] . Therefore, any of these sectors provide separable 
Hilbert spaces contained in Ti^:^^. We denote them as ?^e,Aj,A| = He^T-Lx* ^Tix*. Note that the 
removal of the states with support in f = means that there is no analog of the cosmological 
singularity (classically located in pi = 0) in our quantum theory. We thus solve the singularity 
already at the level of superselection in a very simple way. 

We then restrict the study to any of these superselection sectors, and expand a general 
solution in the form {ip] = Y.yeC+ Sc^.eW, T^Q^eW, 4'{v,^^eK^^e^s){^^^eK^^e>'s\- obtain 
that the constraint (V^jC^ = leads to a recurrence equation that relates the combination of 
states 

/ v + A \ ( + 4 A / f; + 2 

+ 4, A^,A5) = V'U + 4,A^, ^y-p^A5 j +VU + 4, — p^A^,A5j + + 4, — ^A^, A^ 

+ V U + 4,A^, A5 +^ U + 4, A^,^— -A5 +VU + 4, ^— A,,, A^ 

V V I \ V V + 2 J \ V + 2 V 

(3) 

with data in the previous sections v and v — A. Therefore, if we regard v as an internal 
time the constraint can be interpreted as an evolution equation in it. It turns out that, 
when particularized to the initial time v = e, the constraint simply gives the combination 
of states V'+(£ + 4, Ao-jA^) in terms of some initial data ^{e, X'^, X'^). From the combinations 
tp+ie + 4, Ao-, A^), it is possible to determine any of the individual terms ip{e + 4, A^, A^), since 
the system of equations that relate the former ones with the latter ones is formally invertible. 
This issue is very non-trivial. Actually, both the separability of the support of A^ and the fact 
that it is dense in M"*" have been essential to show it (the details can be found in [^). Therefore, 
the initial value problem is well-posed and it is indeed feasible to make the above interpretation 
of the constraint as an evolution equation in v. 

In conclusion, the physical solutions of the Hamiltonian constraint are completely determined 
by the set of initial data {ip{e, X^^, Xs) = '4'{e,uj£X*,Os'^s)' ^ei^s G VVe}, and we can identify 
solutions with this set. We can also characterize the physical Hilbert space as the Hilbert space 



of the initial data. In order to endow the set of initial data with a Hilbert structure, one can take a 
complete set of observables forming a closed algebra, and impose that the quantum counterpart 
of their complex conjugation relations become adjointness relations between operators. This 
determines a unique (up to unitary equivalence) inner product. 

Before doing that, it is most convenient changing the notation. Let us introduce the variables 
Xa = In(Aa) = ln(A*) + pe- Note that pe takes values in a dense set of the real line, given by the 
logarithm of the points in the set Ws- We will denote that set by -Z^. Then, a set of observables 
acting on the initial data ^p{xa,xs) ■= ipi^j Xa,xs) is that formed by the operators e*^" and US" , 
with Pa G and a = a,6, defined as 

e*^-'0(j;^,X5) = e*''"-0(x^,X5), C/^"-0(x^, x^) = + p^, x^), (4) 

and similarly for e"* and Ug^ . These operators provide an overcomplete set of observables and 
are unitary in Tix* (8> 'Hx*, according with their reality conditions. Therefore, we conclude that 
this Hilbert space is precisely the physical Hilbert space of the vacuum Bianchi I model. 

3. Conclusions 

Adopting the improved dynamics of reference [8], we have completed the quantization of the 
Bianchi I model in LQC providing the physical Hilbert space and a complete set of physical 
observables, task so far not achieved. 

With that aim, first we have needed to determine the structure of the superselection sectors 
that arise in the quantum model owing to the polymeric representation of the geometry. These 
sectors display an involved structure concerning the anisotropies, as a consequence of the 
complicated action of the Hamiltonian constraint on those variables, inherent to the adopted 
improved dynamics. Actually, the values of the anisotropies run over a dense set of the positive 
real line that turns out to be countable. In contrast, the volume takes values in simple 
semilattices of constant step. As a result, the superselection sectors are separable. 

The Hamiltonian constraint provides a difference equation in the volume, and can be 
interpreted as an evolution equation, being the volume the variable playing the role of the 
internal time. The separability of the support of the anisotropies and the fact that it is dense in 
R"*" are essential to show that the associated initial value problem is well-posed, and therefore the 
above interpretation is valid. Since the initial data determines the solution, we can characterize 
the physical Hilbert space as the Hilbert space of the initial data. The physical inner product 
can be determined by imposing reality conditions on a complete set of observables. The result 
is that the physical Hilbert space coincides with the tensor product of the superselection sectors 
of the two anisotropy variables. 
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